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DOMAIN  DECOMPOSITION  ALGORITHMS 
WITH  SMALL  OVERLAP 

MAKSYMILIAN  DRYJA  *   AND  OLOF  B.  WIDLUND  ' 

Abstract.  Numerical  experiments  have  shown  that  two-level  Schwarz  methods  often  perform  very 
well  even  if  the  overlap  between  neighboring  subregions  is  quite  small.  This  is  true  to  an  even  greater 
extent  for  a  related  algorithm,  due  to  Barry  Smith,  where  a  Schwarz  algorithm  is  applied  to  the  reduced 
linear  system  of  equations  that  remains  after  that  the  variables  interior  to  the  subregions  have  been 
eliminated.  In  this  paper,  a  supporting  theory  is  developed. 

Key  words,  domain  decomposition,  elliptic  finite  element  problems,  preconditioned  conjugate 
gradients,  Schwarz  methods 
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1.  Introduction.  Over  the  last  decade,  a  considerable  interest  has  developed 
in  Schwarz  methods  and  other  domain  decomposition  methods  for  partial  differential 
equations;  cf.  e.g.  the  proceedings  of  five  international  symposia  [34,17,18,35,19].  A 
general  theory  has  evolved  and  a  substantial  number  of  new  algorithms  has  been  de- 
signed, analyzed  and  tested  numerically.  Among  them  are  two-level,  additive  Schwarz 
methods  first  introduced  in  1987;  cf.  Dryja  and  Widlund  [28,25,29,30,55].  For 
related  work  see  also  BJ0rstad,  Moe,  and  Skogen  [1,2,3],  Cai  [10,11,12],  Mathew 
[40,42,41],  Matsokin  and  Nepomnyaschikh  [43],  Nepomnyaschikh  [44],  Skogen  [47], 
Smith  [48,52,49,50,51],  and  Zhang  [59,60].  As  shown  in  Dryja  and  Widlund  [30],  a 
number  of  other  domain  decomposition  methods,  in  particular  those  of  Bramble,  Pas- 
ciak,  and  Schatz  [5,6],  can  also  be  derived  and  analyzed  using  the  same  framework. 
Recent  efforts  by  Bramble,  Pasciak,  Wang,  and  Xu  [7],  and  Xu  [56]  have  extended 
the  general  framework  making  a  systematic  study  of  multiplicative  Schwarz  meth- 
ods possible.  The  multiplicative  algorithms  are  direct  generalizations  of  the  original 
alternating  method  discovered  more  than  120  years  ago  by  H.A.  Schwarz  [46].  For 
other  current  projects,  which  also  use  the  Schwarz  framework,  see  Dryja,  Smith,  and 
Widlund  [27]  and  Dryja  and  Widlund  [32,33]. 

When  a  two-level  method  is  used,  the  restrictions  of  the  discrete  elliptic  problem 
to  overlapping  subregions,  into  which  the  given  region  has  been  decomposed,  are 
solved  exactly  or  approximately.  These  local  solvers  form  an  important  part  of  a 
preconditioner  for  a  conjugate  gradient  method.  In  addition,  in  order  to  enhance  the 
convergence  rate,  the  preconditioner  includes  a  global  problem  of  relatively  modest 
dimension. 
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Generalizations  to  more  than  two  levels  have  also  been  developed;  see  e.g.  Bram- 
ble, Pasciak,  and  Xu  [8],  Dryja  and  Widlund  [31],  and  Zhang  [58,59,60];  here  the 
families  of  domain  decomposition  methods  and  multigrid  algorithms  merge.  Recently 
there  has  also  been  a  considerable  interest  in  nonsymmetric  and  indefinite  problems; 
cf.  e.g.  Bramble,  Leyk,  and  Pasciak  [4],  Cai  [10,11,12],  Cai,  Gropp,  and  Keyes  [13], 
Cai  eind  Widlund  [14,15],  Cai  and  Xu  [16],  and  Xu  [57].  However,  in  this  paper,  we 
work  exclusively  with  two-level  methods  for  positive  definite,  symmetric  problems. 

The  main  result  of  our  eeurly  study  of  two-level  Schwarz  methods  shows  that  the 
condition  number  of  the  operator,  which  is  relevant  for  the  conjugate  gradient  itera- 
tion, is  uniformly  bounded  if  the  overlap  between  neighboring  subregions  is  sufficiently 
generous  in  proportion  to  the  diameters  of  the  subregions. 

Our  current  work  has  been  inspired  very  directly  by  several  series  of  numerical 
experiments  that  indicate  that  the  rate  of  convergence  is  quite  satisfactory  even  for  a 
small  overlap  and  that  the  running  time  of  the  programs  is  often  the  smallest  when 
the  overlap  is  at  a  minimum.  The  number  of  conjugate  gradient  iterations  is  typically 
higher  in  such  a  case  but  this  can  be  compensated  for  by  the  fact  that  the  local 
problems  are  smaller  and  therefore  cheaper  to  solve;  cf.  in  particular  Bj0rstad,  Moe, 
and  Skogen  [2],  Bj0rstad  and  Skogen  [3],  Cai  [10,11],  Cai,  Gropp,  and  Keyes  [13],  and 
Skogen  [47].  If  the  local  problems  are  themselves  solved  by  an  iterative  method,  then  a 
smaller  overlap  will  give  better  conditioned  local  problems  and  therefore  a  higher  rate 
of  convergence;  see  Skogen  [47]  for  a  detailed  discussion  of  this  effect.  All  this  work 
also  shows  that  these  algorithms  are  relatively  easy  to  implement.  Recent  experiments 
by  Gropp  and  Smith  [37]  for  problems  of  linear  elasticity  provide  strong  evidence  that 
these  methods  can  be  quite  effective  even  for  very  ill-conditioned  problems.  In  this 
paper,  we  show  that  the  condition  number  of  the  preconditioned  operator  for  the 
algorithm,  introduced  in  1987  by  Dryja  and  Widlund  [28],  is  bounded  from  above  by 
const.{l  +  (H/S)).  Here  H  measures  the  diameter  of  a  subregion  and  6  the  overlap 
between  neighboring  subregions.  We  note  that  H/6  is  a  measure  of  the  aspect  ratio 
of  the  subregion  common  to  two  overlapping  neighboring  subregions. 

We  then  turn  our  attention  to  a  very  interesting  method,  introduced  in  1989  by 
Barry  Smith  [52,48].  It  is  known  as  the  vertex  space  (or  Copper  Mountain)  algorithm. 
Numerical  experiments,  for  problems  in  the  plane,  have  shown  that  this  method  con- 
verges quite  rapidly  even  for  problems,  which  were  originally  very  ill-conditioned,  even 
if  the  overlap  is  very  modest;  cf.  Smith  [48].  For  additional  work  on  variants  of  this 
method,  see  Chan  and  Mathew  [20,21],  Chan,  Mathew,  and  Shao  [22]. 

When  Smith's  algorithm  is  used,  the  given  large  linear  system  of  algebraic  equa- 
tions, resulting  from  a  finite  element  discretization  of  an  elliptic  problem,  is  first 
reduced  in  size  by  eliminating  all  variables  associated  with  the  interiors  of  the  nonover- 
lapping  substructures,  {^,},  into  which  the  region  has  been  subdivided.  The  reduced 
problem  is  known  as  the  Schur  complement  system  and  the  remaining  degrees  of  free- 
dom are  associated  with  the  set  {dO.,]  of  substructure  boundaries,  which  form  the 
interface  T  between  the  substructures.  The  preconditioner  of  this  domain  decompo- 
sition method,  classified  as  a  Schwarz  method  on  the  interface  in  Dryja  and  Widlund 


[30],  is  constructed  from  a  coarse  mesh  problem,  with  the  substructures  serving  as 
elements,  and  a  potentially  large  number  of  local  problems.  The  latter  correspond 
to  an  overlapping  covering  of  F,  with  each  subset  corresponding  to  a  set  of  adjacent 
interface  variables. 

Smith's  main  theoretical  result,  given  in  [52,48],  is  quite  similar  to  that  for  the 
original  two-level  Schwarz  method;  the  condition  number  of  this  domain  decompo- 
sition algorithm  is  uniformly  bounded  for  a  class  of  second  order  elliptic  problems 
provided  that  there  is  a  relatively  generous  overlap  between  neighboring  subregions 
that  define  the  subdivision  of  the  domain  decomposition  method.  In  this  paper,  we 
show  that  the  condition  number  of  the  iteration  operator  grows  only  in  proportion  to 
(1  +  \og{H/6)y.  We  note  that  even  for  a  minimal  overlap  of  just  one  mesh  width  /i, 
this  bound  is  as  strong  as  those  for  the  well  known  iterative  substructuring  methods 
considered  by  Bramble,  Pasciak,  and  Schatz  [5,6],  Dryja  [24],  Dryja,  Proskurowski, 
and  Widlund  [26],  Smith  [49],  and  Widlund  [54];  cf.  also  Dryja,  Smith,  and  Widlund 
[27].  We  also  note  that  the  successful  iterative  substructuring  methods  for  problems 
in  three  dimensions,  require  the  use  of  more  complicated  coarse  subspaces  and  that 
therefore  Smith's  method,  considered  in  this  paper,  seems  to  offer  axi  advantage. 

2.  Some  Schwarz  Methods  for  Finite  Element  Problems.  As  usual,  we 
write  our  continuous  and  finite  element  elliptic  problems  as:  Find  u  £  V,  such  that 

a{u,v)  =  f{v),  W  V  e    V, 

and,  find  u^  E  V",  such  that 

(1)  a{u,,v,)  =  f{v,),'^v,e   V\ 

respectively.  We  assume  that  the  bilinear  form  a{u,v)  is  selfadjoint  and  elliptic  and 
that  it  is  bounded  in  F  x  K  In  the  case  of  Poisson's  equation,  the  bilinear  form  is 
defined  by 

(2)  a{u,v)  =  f  Vu-Vvdx  . 

Jn 

We  cissume  that  fi  is  a  Lipschitz  region  in  i?",n  =  2,3,  and  that  its  diameter  is  on 
the  order  of  1.  (We  will  follow  Neceis  [45]  when  defining  Lipschitz  regions  and  Sobolev 
spaces  on  Q,.)  The  bilinear  form  a(u,  v)  is  directly  related  to  the  Sobolev  space  H^{^) 
that  is  defined  by  the  semi-norm  and  norm 

l"l//>(n)  =°("'^)   a^^   M\m(Q)  =  W\HHQ)  +  h\\h{Q)^ 

respectively.  When  we  axe  considering  subregions  fi,  C  Cl,  of  diameter  H,  we  use  a 
different  relative  weight,  obtained  by  dilation, 

|2 


ll"llH'(n.)  =  l"lH'(n.)  +  ;^ll"lli:^(n.)- 


Whenever  appropriate,  we  tacitly  assume  that  the  elements  of  ifo(ri,),  the  subspace 
of  H^{Ct,)  with  zero  trace  on  the  boundary  dQ,,  are  extended  by  zero  to  Q  \  f2,. 
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To  avoid  unnecessary  complications,  we  confine  our  discussion  to  Poisson's  equa- 
tion, to  homogeneous  Dirichlet  conditions  and  to  continuous,  piecewise  linear  finite  ele- 
ments and  a  polygonal  region  Q.  It  is  well  known  that  the  resulting  space  V'^  C  Hq{Q,), 
i.e.  it  is  conforming. 

For  the  problem  considered  in  a  1870  paper  by  H.A.  Schwarz  [46],  two  overlapping 
subregions,  Q.[  and  fij  axe  used;  the  union  of  the  two  is  Q,.  There  are  two  sequential, 
fractional  steps  of  the  iteration  in  which  the  approximate  solution  of  the  elliptic  equa- 
tion on  Q,  is  updated  by  solving  the  given  problem  restricted  to  the  subregions,  one 
at  a  time.  The  most  recent  values  of  the  solution  are  used  as  boundary  values  on  the 
part  of  dCt'-,  the  boundary  of  Q.[  that  is  not  a  part  of  dQ. 

The  finite  element  version  of  the  algorithm  can  conveniently  be  described  in  terms 
of  projections  P,  :  V^  -^  V/"  =  H^(n',)f]V\  defined  by 

(3)  a{PiVh,  <t>h)  =  aivH,  (f>h),  V.^,  G  V;'' . 

It  is  easy  to  show  that  the  error  propagation  operator  of  this  multiplicative  Schwarz 
method  is 

[I  -  P,){I  -  P,). 
This  algorithm  can  therefore  be  viewed  as  a  simple  iterative  method  for  solving 

{Pi  +  P2  -  P2Pi)uh  =  gn, 

with  an  appropriate  right-hand  side  g/,. 

This  operator  is  a  polynomial  of  degree  two  and  therefore  not  ideal  for  parallel 
computing  since  two  sequential  steps  are  involved.  This  effect  is  further  pronounced 
if  more  than  two  subspaces  are  used.  Therefore,  it  is  often  advantageous  to  collect 
subregions,  which  do  not  intersect,  into  groups;  the  subspaces  of  each  group  can  then 
be  regarded  as  one.  The  number  of  subspaces  is  thus  reduced,  and  the  algorithm 
becomes  easier  to  parallelize.  Numerical  experiments  with  multiplicative  Schwarz 
methods  have  also  shown  that  the  convergence  rate  often  is  enhanced  if  such  a  strat- 
egy is  pursued;  this  aproach  is  similar  to  a  red-black  or  multi-color  ordering  in  the 
context  of  classical  iterative  methods.  In  the  case  of  an  additive  Schwarz  methods, 
this  ordering  only  serves  as  a  device  to  facilitate  the  analysis. 

In  the  additive  form  of  the  algorithm,  we  work  with  the  simplest  possible  poly- 
nomial of  the  projections:  The  equation 

(4)  Puk  =  {Pi+P2  +  ---  +  PN)uk=g[, 

is  solved  by  an  iterative  method.  Here  P,  :V  —^  V,,  and  V  =  Vi -\ hV^.  Since  the 

operator  P  can  be  shown  to  be  positive  definite,  symmetric,  with  respect  to  a(-,  •),  the 
iterative  method  of  choice  is  the  conjugate  gradient  method.  Equation  (4)  must  also 
have  the  same  solution  as  equation  (1),  i.e.  the  correct  right-hand  side  must  be  found. 
This  can  easily  be  arranged;  see  e.g.  [29,30,55].  Much  of  the  work,  in  particular  that 
which  involves  the  individual  projections,  can  be  carried  out  in  parallel. 
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2.1.  The  Dryja-Widlund  Algorithm.  We  now  describe  the  special  additive 
Schwarz  method  introduced  in  Dryja  and  Widlund  [28];  cf.  also  Dryja  [25]  and  Dryja 
and  Widlund  [29].  We  start  by  introducing  two  triangulations  of  Q.  into  nonover- 
lapping  substructures  fi,  and  into  elements.  We  obtain  the  elements  by  subdividing 
the  substructures.  We  cdways  assume  that  the  two  triangulations  are  shape  regular, 
cf.  e.g.  Ciarlet  [23],  and,  to  simpHfy  our  arguments,  that  the  diameters  of  all  the 
substructures  are  on  the  order  of  H.  In  this  algorithm,  we  use  overlapping  subregions 
obtained  by  extending  each  substructure  Q.,  to  a  larger  region  fi(.  The  overlap  is  said 
to  be  generous  if  the  distance  between  the  boundaries  dQ,  and  dQ,[  is  bounded  from 
below  by  a  fixed  fraction  of  H.  We  always  assume  that  dCl',  does  not  cut  through  any 
element.  We  carry  out  the  same  construction  for  the  substructures  that  meet  the 
boundary  except  that  we  cut  off  the  part  of  fi^  that  is  outside  of  f2. 

We  remark  that  other  decompositions  are  also  of  interest.  In  particular,  the 
analysis  in  Section  4  extends  immediately  to  the  case  when  no  degrees  of  freedom  axe 
shared  between  neighboring  subregions.  In  this  case  the  distance  between  dQ.[  and 
dQ'j  is  just  h  for  neighboring  subregions.  This  additive  Schwcirz  method  corresponds 
to  a  block  Jacobi  preconditioner  augmented  by  a  coarse  solver. 

For  this  Schwarz  method,  the  finite  element  space  is  represented  as  the  sum  of 
N+1  subspaces 

y''  =  Fj>  +  v;''  +  •  •  •  +  V;^ . 

The  first  subspace  V^''  is  equal  to  V^ ,  the  space  of  continuous,  piecewise  linear  func- 
tions on  the  coarse  mesh  defined  by  the  substructures  Q,.  The  other  subspaces  are 
related  to  the  subdomains,  in  the  same  way  as  in  the  original  Schwarz  algorithm,  i.e. 

It  is  often  more  economical  to  use  approximate  rather  than  exact  solvers  for  the 
problems  on  the  subspaces.  The  approximate  solvers  are  described  in  the  following 
terms:  Let  b,{u,v)  be  an  inner  product  defined  on  V^'^  x  V^''  and  assume  that  there 
exists  a  constant  u  such  that 

(5)  a(u,u)  <  u;fe.(u,u)  ,Vu  e  Vf  . 

In  terms  of  matrices,  this  inequality  becomes  a  one-sided  bound  of  the  stiffness  matrix, 
corresponding  to  a(-,  •)  and  V^'',  in  terms  of  the  matrix  corresponding  to  the  bilinear 
form  6,(u,  u). 

An  operator  T,  :  V'^  — >  V/^,  which  replaces  P,,  is  now  defined  by 

(6)  6.(T,u,,^0  =  a(",<^/.),  V0,6V;.^ 

It  is  easy  to  show  that  the  operator  T,  is  positive  semidefinite  and  symmetric  with 
respect  to  a(-,-)  and  that  the  minimal  constant  u.'  in  equation  (5)  is  UTiHa-  Additive 
and  multiplicative  Schwarz  methods  can  now  be  defined  straightforwardly  in  terms  of 
polynomials  of  the  operators  T,.  We  note  that  if  exact  solvers,  and  thus  the  projections 
P,,  are  used,  then  u;  =  1. 
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2.2.  Smith's  Algorithm.  Smith's  method  has  previously  been  described  in 
Smith  [52,48].  Let  K  be  the  stiffness  matrix  given  by  the  bihnear  form  of  (2).  In 
the  first  step  of  this,  and  many  other  domain  decomposition  methods,  the  unknowns 
of  the  Hneax  system  of  equations 

Kx  =  6, 

which  correspond  to  the  the  interiors  of  the  substructures  are  ehminated.  We  now 
describe  this  procedure  in  some  detail. 

Let  i\  <')  be  the  stiffness  matrix  corresponding  to  the  bihnear  auXuh^Vh)  which  rep- 
resents the  contribution  of  the  substructure  fi,  to  the  integral  aQ{uh,Vh)  =  a{uhiVh)- 
Let  X  cind  y  be  the  vectors  of  nodal  values  that  correspond  to  the  finite  element  func- 
tions Uh  and  Vk,  respectively.  Then  the  stiffness  matrix  K  of  the  entire  problem  can 
be  obtained  by  using  the  method  of  subassembly  defined  by  the  formula 

Here  x^'^  is  the  subvector  of  nodal  parameters  associated  with  f2,,  the  closure  of  f2,. 
We  represent  A'^''  as 

(   A'l?     a!2  \ 

Here  we  have  divided  the  subvector  x*''  into  two,  x/  and  x^  ,  corresponding  to  the 
variables  which  are  interior  to  the  substructure  and  those  which  are  shared  with  other 
substructures,  i.e.  those  associated  with  the  nodal  points  of  5^2, .  Since  the  interior 
variables  of  Q,  are  coupled  only  to  other  variables  of  the  same  substructure,  they 
can  be  eliminated  locally  and  in  parallel.  The  resulting  reduced  matrix  is  a  Schur 
complement  and  is  of  the  form 

(7)  5C)  =  Aa-Ai2^Aj?-Al2. 

From  this  follows  that  the  Schur  complement,  corresponding  to  the  global  stiffness 
matrix  A',  is  given  by  S  where 

(8)  xlSys  =  E-B^S^"yB- 

If  the  local  problems  are  solved  exactly,  what  remains  is  to  find  a  sufficiently  accurate 
approximation  of  the  solution  of  the  linear  system 

(9)  SXB  =  9B. 

It  is  convenient  to  rewrite  (9)  in  variational  form.  Let  s,{uh,Vh)  and  s{uh,Vh} 
denote  the  forms  defined  by  (7)  and  (8),  respectively,  i.e. 

s,iuh,Vh)  =  x^'S<'>y^'       and      s{uh,Vh)  =  x^SyB- 
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Equation  (9)  can  then  be  rewritten  as 

(10)  s{uh,VH)  =  ig,VH)mr),  yvh  e  v\r). 

EeTeT  =  \Jdn,\dQ. 

Problem  (10)  will  be  solved  by  an  iterative  method  of  additive  Schwaxz  type.  The 
most  importemt  difference  between  this  algorithm  and  that  of  the  previous  subsection 
is  that  we  are  now  working  with  the  trace  space  H^I'^{T)  instead  ol  H^{Vi);  see  Section 
3  for  a  definition  of  H^I'^{T). 

It  is  well  known  that 

(11)  x^)^5<'\^>  =  minx(')^A'(')x('). 

Therefore,  if  uh  is  the  minimsd,  the  discrete  harmonic  extension  of  the  boundary  data 
represented  by  ib,  then 

Smith's  algorithm  can  now  be  described  in  terms  of  a  subspace  decomposition. 
We  use  the  same  coarse  space  as  in  the  previous  subsection,  i.e.  V^ ,  but  we  restrict 
its  values  to  T.  In  the  case  when  the  originad  problem  is  two  dimensional,  we  introduce 
one  subspace  for  each  interior  edge  and  one  for  each  vertex  of  the  substructures.  An 
edge  space  is  defined  by  setting  all  nodal  values,  except  those  associated  with  the 
interior  of  the  edge  in  question,  to  zero  .  Similarly,  a  vertex  space  is  obtained  by 
setting  to  zero  all  values  at  the  nodes  on  F,  which  are  at  a  distance  greater  than  6. 
For  many  more  details  and  a  discussion  of  implementation  details,  see  Smith  [52,48]. 

In  the  case  when  the  original  problem  is  three  dimensional,  we  introduce  one 
subspace  for  each  interior  face,  edge,  and  vertex.  The  elements  of  a  face  subspace 
vanish  at  cdl  nodes  on  F  that  do  not  belong  to  the  interior  of  the  face.  Similarly,  an 
edge  space  is  supported  in  the  strips  of  width  6,  which  belong  to  the  faces,  which  have 
this  edge  in  common.  Finally,  a  vertex  space  is  defined  in  terms  of  the  nodes  on  F 
that  are  within  a  distance  S  of  the  vertex. 

2.3.  Basic  Theory.  In  order  to  estimate  the  rate  of  convergence  of  our  special, 
or  any  other,  additive  Schwaxz  methods,  we  need  upper  and  lower  bounds  for  the 
spectrum  of  the  operator  relevant  in  the  conjugate  gradient  iteration.  A  lower  bound 
can  be  obtziined  by  using  the  following  lemma;  cf.  e.g.  Dryja  and  Widlund  [29,33]  or 
Zhang  [58]. 

Lemma  1.  Let  Ti  be  the  operators  defined  in  equation  (6)  and  let  T  =  To  +  Ti  + 
•  •  •  +  Tf4.  Then 

a{T~^u,u)  =    min   ^6,(u,,u,),  u,  G  V,. 

Therefore,  if  a  representation,  u  =  Y^Ui,  can  be  found,  such  that 

Y,b,{u,,u,)<CMu,u),    VueV\ 


then 

An  upper  bound  for  the  spectrum  of  T  is  often  obtained  in  terms  of  strength- 
ened Cauchy-Schwarz  inequeilities  between  the  different  subspaces.  Note  that  we  now 
exclude  the  index  0;  the  coarse  subspace  is  treated  separately. 

Definition  1.  The  matrix  £  =  {e,;}  is  the  matrix  of  strengthened  Cauchy- 
Schwarz  constants,  i.e.  dj  is  the  smallest  constant  for  which 

(12)  Hv„v,)\  <  £o|l^.||a||^,l|a  yv,  €  V,  yv,  G  V,  ,i,j  >  1  , 

holds. 

The  following  lemma  is  easy  to  prove;  cf.  Dryja  and  Widlund  [33]. 

Lemma  2.  Let  p{S)  he  the  spectral  radius  of  the  matrix  S.  Then,  the  operator  T 
satisfies, 

T<u:ip{e)  +  l)I. 

For  the  pcirticular  algorithms  considered  in  this  paper,  it  is  very  easy  to  show  that 
there  is  a  uniform  upper  bound.  In  fact  by  collecting  and  merging  local  subspaces 
that  belongs  to  nonoverlapping  subregions,  the  number  of  subspaces,  and  piS),  can 
be  made  uniformly  bounded;  see  Section  4  for  an  alternative  argument. 

By  combining  Lemma  1  and  2,  we  obtain 

Theorem  1.  The  condition  number  k(T)  of  the  operator  T  of  the  additive 
Schwarz  method  satisfies, 

k{T)  =  Xr.a.{T)/\m,niT)  <  u;(/>(^)  +  l)C^    . 

In  the  multiplicative  case,  we  need  to  provide  an  upper  bound  for  the  spectral 
radius,  or  norm,  of  the  error  propagation  operator 

(13)  Ej  =  iI-Tj)---{I-To). 

The  following  theorem  is  a  variant  of  a  result  of  Bramble,  Pasciak,  Wang  and  Xu  [7]; 
cf.    also  Cai  and  Widlund  [15],  Xu  [56]  or  Zhang  [58].    Note  that  this  bound  is  also 
given  in  terms  of  the  same  three  pEirameters  that  appears  in  Theorem  1 . 
Theorem  2.  In  the  symmetric,  positive  definite  case, 


'^J\\a   < 


1_  (2-a;) 


M         (2wV(^)^  +  l)C2 


We  note  that  this  formula  is  useless  if  u;  >  2;  since  ||7  -  T,\\a  >  1  if  ||r,||a  >  2,  the 
assumption  that  u;  <  2  is  most  natural.  If  we  wish  to  use  a  multiplicative  algorithm 
ajid  u  is  too  large,  we  can  sccde  the  bilinear  forms  fe,(-,  •)  suitably. 

In  this  paper,  our  results  Eire  only  formulated  for  additive  algorithms.  The  corre- 
sponding bounds  for  the  multiplicative  variants  can  easily  be  worked  out  as  applica- 
tions of  Theorem  2. 
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3.  Technical  Tools.  In  this  section,  we  collect  a  number  of  techniced  tools  that 
are  used  to  prove  our  main  results.  Some  of  these  tools  axe  quite  familiar  to  specialists 
of  the  field.  Others  have,  to  our  knowledge,  not  previously  been  used  in  the  analysis 
of  domain  decomposition;  cf.  Il'in  [38]  for  some  similar  inequalities. 

As  before,  J)  C  /2",  n  =2  or  3,  is  a  bounded,  polygonal  region  and  {f2,}  a  nonover- 
lapping  decomposition  of  Q,  into  substructures.  To  simplify  our  considerations,  we  now 
assume  that  the  substructures  are  squares  or  a  cubes;  cf.  e.g.  Necas  [45]  where  sim- 
ple maps  and  paxtions  of  unity  axe  used  to  derive  bounds  for  Lipschitz  regions  from 
bounds  for  such  special  regions;  if  we  can  handle  a  corner  of  a  square  or  cube,  then  we 
can  analyze  the  general  polygonal  case.  Our  estimates,  given  in  the  next  two  sections, 
are  developed  for  one  substructure  at  time;  our  arguments  can  be  modified  to  make 
them  valid  for  any  shape  regular  substructure  with  a  boundary  consisting  of  a  finite 
number  of  smooth  curves  or  surfaces. 

As  before,  let  T  =  \JdQ,i  \  d^.  Let  F^,,  C  ^i  be  the  set  of  points  that  is  within  a 
distance  6  of  F. 

Lemma  3.  Lei  u  be  an  arbitrary  element  of  H^{Q,i).  Then, 

ll"lli=(r,.)  <  C8\{1  +  ^/6)|u|^,,„.)  +  l/^<5||u||i.(„,,)  . 

Proof.  We  first  consider  a  square  region  (0,  H)  x  (0,  H)  in  detail;  the  extension  of 
the  proof  to  the  Ccise  of  three  dimensions  is  straightforward.  Since 

u{x, 0)  =  u(x, y)  -  J    — Q^ — dr  , 

we  find,  by  elementary  axguments,  that 

H        \u{x,0)\'dx<2        /     \u{x,y)\'dxdy  +  H'        /     \^\'dxdy  . 
Jo  Jo   Jo  Jo   Jo      oy 

Therefore, 


H  j"  \u{xM^dx  <  2||K||i,(^.)  +  i?>|2,, 


(n.) 


Now  consider  the  integral  over  a  narrow  subregion  next  to  one  of  the  sides  of  the 
square.  Using  similar  arguments,  we  obtain 

(14)  j"  f^   \u{x,y)\^dxdy  <  S'\u\jf^^^^^  +  26  j^  |u(x,0)|2(ix  . 

By  combining  this  and  the  previous  inequality,  we  obtain 

rH      rS  0 

Jo    Jo  '"(^'^)l''^^^^-'^'l"lw>(n.)  +  2<5(-^l|i.||i,(n.)  +  i/|upH>(f^,))' 

as  required. 

The  modifications  necessary  for  the  case  of  an  arbitrary,  shape  regulax  substruc- 
ture and  the  extension  of  the  proof  to  the  case  of  three  dimensions  are  routine.    D 
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Lemma  4.  Let  Uh  be  a  continuous,  piecewise  quadratic  function  defined  on  the 
finite  element  triangulation  and  let  IhUh  G  V^  be  its  piecewise  linear  interpolant  on 
the  same  mesh.    Then  there  exists  a  constant  C,  independent  of  h  and  H,  such  that 

UhUhlnHQ.)  <  C\uh\HHu,)  • 

The  same  type  of  bounds  hold  for  the  L2iH^'^,  and  Hqq     norms  and  it  can  also  be 
extended,  with  different  constants,  to  the  case  of  piecewise  cubic  functions,  etc. 
Proof.  It  is  elementary  to  show  that, 

\huh\H^(Q,)  <  2{\IhUh  -  UhlnnQ,)  +  \^h\HHQ,))- 

Consider  the  contribution  to  the  first  term  on  the  right  hand  side  from  an  individual 
element  K.  We  obtain 

\IhUh  -  UklnHK)  ^  Ch^\uh\HHf^')  <  C\uh\Hi{K) 

by  using  a  standard  error  bound  and  an  elementary  inverse  inequality  for  quadratic 
polynomials.  The  bound  in  L2  follows  from  the  linear  independence  of  the  standard 
finite  element  basis  for  the  space  of  quadratic  polynomials,  cf.  Ciarlet  [23].  The 
bounds  for  the  other  norms,  which  are  defined  below,  can  be  obtained  by  interpolation 
in  Sobolev  spaces;  cf.  e.g.  Lions  and  Magenes  [39].    D 

We  now  turn  to  the  other  auxiliary  results  needed  in  the  analysis  of  Smith's 
algorithm.  We  begin  by  providing  an  expression  for  the  norm  of  H^^^;  cf.  Chapters 
1.3.2  and  1.5  of  Grisvard  [36]  for  a  detailed  discussion. 

Let  I  Q  E}  be  an  open  interval  of  diameter  H.  Then, 


(15)  ||„|i;,.„,„  =  jf^Mfi^,,,,  +  ^ll„l 


L2{iy 


The  relative  weight  of  the  two  terms  is  obtained,  by  dilation,  from  the  norm  defined 

on  a  region  of  diameter  1. 

It  is  well  known  that  the  extension  by  zero  of  the  elements  of  H^I^{I)  does  not 

define  a  continuous  map  into //''/^(i?^);  cf.  Lemma  1.3.2.6  of  Grisvard  [36]  or  Lions  and 

Magenes  [39].  The  largest  subspace  for  which  this  extension  operator  is  continuous  is 

1/2 
-^00  (■^)i  which  is  defined  in  terms  of  the  norm  obtained  by  replacing  the  last  term  of 

(15)  by 


(16)  / 


d[s) 


ds. 


Here  d{s)  is  the  distance  to  the  end  points  of  /. 

In  the  case  of  a  subset  ^  of  the  boundary  of  a  three  dimensional  region,  the 
formula  (15)  is  valid  after  replacing  \s  —  <|^  by  \s  —  t\^  and  /  by  '^ .  However,  for 
our  purposes,  it  is  more  convenient  to  use  an  alternative  formula;  cf.   Dryja  [24]  and 
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Lemma  5.3,  Chapter  2  of  Necas  [45].  In  the  special  case  of  a  square  with  side  H,  the 
semi-norm  is  defined  by 

^     '     Jo    Jo  \s\-U\^  Jo    Jo  \s2-t2V 

To  obtain  the  norm  for  the  subspace  Hll^{{0,Hf),  we  add  a  weighted  norm 

Jo    Jo  d{s) 


just  as  in  (16). 

In  addition  to  the  space  F'',  we  will  also  use  a  coarser  space  V^,  defined  on  a 
mesh  with  mesh  size  6,  in  our  proofs.  We  now  formulate  results  that  have  been  used 
extensively  in  work  of  this  kind;  cf.  Dryja  [24],  or  Bramble  and  Xu  [9].  The  first 
inequality  of  the  lemma  is  given  as  Lemma  1  in  Dryja  [24].  The  second  is  part  of  the 
proof  of  his  Lemma  4  in  the  same  paper. 

Lemma  5.  Lei  I  be  an  interval  of  length  H.  Then, 

lh5|lL(/)  <  C-Cl  +  ^og{H/6))\\us\\lu^jy    Vu5  e  v. 
Let  I  be  an  edge  of  a  face  ^  of  diameter  H  of  a  cube.   Then, 

The  next  result  gives  a  bound  which  is  similar  to  the  second  formula  of  Lemma  5. 
However,  the  bound  holds  for  all  of  H^^^. 

Lemma  6.  Let  <if  =  {0,Hy  and  lei  ^s  =  {0,H)  x  (0,6).  Then, 

II"IIL(*,)  <  C6{1  +  log(i//<5))||u||^,/.,^),    Vu  6  H'/\n 

The  same  result  holds  if  we  replace  ^  and  ^^  by  {0,H)  and  (0,6),  respectively. 

Proof.  We  only  provide  a  proof  for  the  first  of  the  two  cases;  the  proof  in  the  other 
case  is  completely  analogous.  Let  Qs  :  H'^{'^)  —^  V\he  the  L2-  projection.  It  is  well 
known  that  ||Q5||//i(*)  is  bounded;  cf.  e.g.  Bramble  and  Xu  [9].  Since,  trivially,  this 
operator  is  also  bounded  in  L2,  it  follows  that  ||<56||//i/2(*)  is  bounded.  By  a  standard 
argument, 

(19)  h  -  Qsu\\l^<t)  <  sc\u\l,„^^y 

We  now  only  need  to  show  that 
(20)  \\Qsu\\l^^^)  <  C6{1  +  \og{H/6mQsn\\lu^^y 

To  prove  (20),  we  use  the  bound  (14)  derived  in  the  proof  of  Lemma  3.  Thus, 

II<3«^^IIL(*,)  <  ^^\Q6u\jf.(^)+26  f     \Qsu{x,0)\^dx  . 
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By  using  an  inverse  inequality,  the  first  term  can  be  replaced  by  ^|Qi"|^i/2(,p)-  The 
second  term  is  estimated  using  Lemma  5.    D 

The  final  lemma  will  be  used  to  estimate  the  weighted  L2  term  in  the  Hqq    norm. 

Lemma  7.  Let  u  e  H^^'^{Q,H).  Then  there  exists  a  constant  C,  such  that 

f"  Mfll!^,  <  c(l  +  \ogiH/6)r\\u\\],UHo,Hy 

•Jo  "S 

Similarly,  let  u  £  H^^^{{0,H)^).  Then  there  exists  a  constant  C,  such  that 

r  ( r  M^i^d.)cf<  <  c(i + iog(i?/<5))^iiHr„,/.((o,^).). 

Jo      Js  s 

Proof.  We  only  consider  the  first  case  in  detail.  Let  Qs  :  H^^^i^)  -^  V*(fi), 
be  the  L2 -projection  onto  the  finite  element  space  with  mesh  size  8.  We  write  u  = 
(u  -  Qsu)  +  Qsu  and  estimate  each  term  separately.  We  first  note  that  by  a  standard 
estimate 

\\u-Qsu\\l<CS\u\l.,,. 

The  bound  for  the  first  term  is  therefore  obtained  by  noting  that  s  >  6  over  the 
interval  of  integration. 

The  other  term  can  be  estimated  by  using  the  first  bound  of  Lemma  5,  which 
results  in  one  logarithmic  factor,  and  the  observation  that, 

from  which  the  second  logarithmic  factor  arises.    D 

4.  Analysis  of  the  Dryja-Widlund  Algorithm.  We  now  use  the  set  F^,,, 
previously  introduced,  to  characterize  the  extent  of  the  overlap.  We  assume  that  all 
X  €  fit,  which  belong  to  at  least  one  additional  overlapping  subregion  0,'^,  lie  in  F^,,-. 

Theorem  3.  In  the  case  when  exact  solvers  are  used  for  the  subproblems,  the 
condition  number  of  the  additive  Schwarz  method  satisfies 

k{P)  <  C(l  +  H/6). 

The  constant  is  independent  of  the  parameters  H,  h  and  S. 

We  note  that  for  the  case  of  two  subregions,  it  is  easy  to  show  that  this  result 
is  sharp.  It  is  routine  to  modify  Theorem  3  to  cover  cases  where  inexact  solvers  axe 
used. 

Proof.  The  proof  is  a  refinement  of  a  result  first  given  in  Dryja  and  Widlund 
[28];  cf.  [29]  for  a  better  discussion.  The  proof  is  equally  valid  for  two  and  three 
dimensions.  We  first  show  that  a  constant  upper  bound  for  the  spectrum  of  P  can  be 
obtained  without  the  use  of  Lemma  2. 

We  note  that  Pi  is  also  an  orthogonal  projection  o{ H^{Q',)f]  V^  onto  Vi.  Therefore, 

a{P,Uh,Uh)  <  aQ'{uh,Uh). 
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Since,  by  construction,  there  is  an  upper  bound,  A'^ci  on  the  number  of  subregions  to 
which  any  z  6  12  can  belong,  we  have 

N 

]^ an; («/,,«/,)  <  Nca{uh,Uh). 
«=i 

In  addition,  we  use  the  fact  that  the  norm  of  Pq  is  equal  to  one  and  obtain 

A„„(P)<(iV,  +  l). 

The  lower  bound  is  obtained  by  using  Lemma  1.  A  natural  choice  of  uq  is  the  La- 
projection  Qn'^h  of  Uh  onto  V^ .  As  previously  pointed  out,  this  projection  is  bounded 
in  Li  as  well  as  H^  and  there  exists  a  constant,  independent  of  h  and  H,  such  that 

(21)  \\u,-QHnh\\L,<CH\\u,\\,. 

Let  Wh  =  Uh  —  QnUh  and  let  u,  =  h{0,Wh)  ,i  =  1,  •  •  • ,  A''.  Here  //,  is  the  inter- 
polation operator  onto  the  space  V^  and  the  6i{x)  define  a  partition  of  unity,  i.e. 
X),  ^i(a:)  =  1.  These  functions  are  chosen  as  nonnegative  elements  of  V^.  It  is  easy  to 
see  that 


Uh 


=  Wh  +^u,. 


In  the  interior  part  of  fi,,  which  does  not  belong  to  Ts,,,  0,  =  I.  This  function 
must  decrecise  to  0  over  a  distance  on  the  order  of  6.  It  is  easy  to  construct  a  partition 
of  unity  with  0  <  ^,  <  1  and  such  that 

IVO.I  <  ^  . 

In  order  to  use  Lemma  1,  we  first  estimate  a{u,,Ui)  in  terms  of  a{wh,tvh)-  We 
consider  the  contribution  from  one  substructure  at  a  time  and  note  that,  trivially, 

Let  K  be  an  element  in  F^,,.  Then,  using  the  definition  of  u,, 

aA-(u.,u.)  <  2aK{0,WH,e,Wh)  +  2aK{hm  "  ^.>'0,  A((^.  "  ^.)u'0)> 

where  ^,  is  the  average  of  ^,  over  the  element  A'.  Using  the  fact  that  the  diameter  of 
A'  is  on  the  order  of  h  and  the  bound  on  V^,,  we  obtain,  after  adding  over  all  the 
relevant  elements, 

ar,..(".,w.)  <  2an.(u)/,,u»/.)  +  -pWh\\mTt,,y 

We  also  need  to  estimate  ar^.Cuj, Uj)  for  the  j  that  correspond  to  neighboring  sub- 
structures. This  presents  no  new  difficulties. 

To  complete  the  proof,  we  need  to  estimate  ||ti'A||x,2(r  )•  We  note  that  each  x  6  12 
is  covered  only  a  finite  number  of  times  by  the  subregions.  We  apply  Lemma  3  to 
the  function  if/,,  sum  over  i  and  use  inequality  (21)  to  complete  the  estimate  of  the 
parameter  Cq  of  Lemma  1 .    D 
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5.  Analysis  of  Smith's  Method.  A  description  of  the  reduction  of  the  original 
Hnear  system  to  one  for  the  degrees  of  freedom  on  F,  and  of  the  algorithm,  have  been 
given  in  Section  2.  We  will  now  work  in  the  H^I'^{T)  norm.  The  fact  that  this  is 
a  weaJcer  norm  than  H^  is  reflected  in  a  stronger  bound  than  that  of  the  previous 
section;  better  bounds  of  the  components  in  the  different  subspaces  can  be  obtained. 
There  are  many  similarities  between  the  two  cases.  Much  of  the  analysis  is  again 
carried  out  one  substructure  at  a  time. 

To  show  that  we  can  work  with  |u/i|^i/2(n  )  instead  of  Xq  5*''x^  ,  we  must  show 
that  these  norms  are  equivalent.  We  use  equation  (11)  and  the  standard  trace  theorem 
to  bound  |u/i|^i/2(Q  x  from  above  by  i^^  5^'^a:g  .  The  proof  of  the  reverse  inequality 
requires  an  extension  theorem  for  finite  element  spaces  given  in  Widlund  [53];  see 
further  discussion  in  Smith  [52]. 

Theorem  4.  In  the  case  when  exact  solvers  are  used  for  the  subproblems,  the 
condition  number  of  the  vertex  space  method  satisfies 

K(P)<C{l  +  \og{H/6)f. 

The  constant  is  independent  of  the  parameters  H,  h  and  8. 

Proof.  As  in  the  proof  of  Theorem  3,  there  is  no  difficulty  in  establishing  a  uniform 
upper  bound  on  the  spectrum  of  P. 

We  now  turn  to  the  lower  bound  in  the  case  where  the  original  problem  is  two 
dimensional  and  thus  the  interface  is  of  dimension  one.  In  order  to  use  Lemma  1,  we 
have  to  decompose  functions  defined  on  T. 

We  use  the  L2(^)  projection  onto  V^  of  the  discrete  harmonic  function  Uk-,  in- 
troduced in  Subsection  2.2,  to  define  the  component  of  the  coarse  space.  We  only 
use  the  values  on  T.  In  addition,  we  use  a  partition  of  unity  to  represent  the  local 
space  components.  In  the  study  of  the  local  spaces,  it  is  sufficient  to  consider  one 
substructure  fi,  at  a  time.  The  partition  of  unity  is  based  on  simple,  piecewise  linear 
functions.  Let  0  <  t  <  H  represent  one  of  the  edges  of  the  boundary  of  this  sub- 
structure and  let  Be{t)  be  a  piece- wise  linear  function,  which  vanishes  for  t  outside 
(0,  //),  grows  linearly  to  1  at  f  =  <5,  is  equal  to  1  for  6  <  t  <  H  —  6  and  drops  to  zero 
linearly  over  the  interval  {H  —  6,  H).  In  the  decomposition,  we  choose  Ih{9eWh)  as  the 
component  corresponding  to  this  edge.  As  in  the  previous  section,  uq  =  Qn^h  and 
u^h  =  Uh  —  Qhuh  is  the  error  of  the  £2—  projection.  It  follows  from  Lemma  4  that  it  is 
sufficient  to  estimate  ||tie||„i/2,n  „,,  where  Ve{t)  =  6At)wk{t).  We  note  that  we  cannot 

use  the  weaker  norm  of  H^^^{0,H)  here;  we  must  estimate  the  H^^'^(dQ,)  norm  of  v^ 
extended  by  zero  to  the  rest  of  the  boundary,  i.e.  ||ue||i,i/2,„  „,. 
We  first  consider 

Jo     Jo  \s  —  t]'' 

and  then  the  additional  term  (16),  which  completes  the  definition  of  the  relevant 
norm.  We  divide  the  interval  [0,H]  into  three  parts,  [0,6],  {6,H  -  6]  and  [H  -  6,H], 
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and  taJce  the  tensor  product  of  [0,  H]  with  itself.  The  double  integral  (22)  is  then  split 
into  the  sum  of  nine.  By  symmetry,  only  six  different  cases  need  to  be  considered. 
The  integral  over  [6,  H  —  S]  x  [S,  H  —  6]  is  completely  harmless.  We  now  consider  the 
diagonal  term  corresponding  to  the  set  [0,<5]  x  [0,^]  and  use  the  identity 

SWH{s)-tWh{t)        {s  +  t){whis)-Whit))    ,    {s-t)(whis)  +  Wh{t)) 
v,{s)  -  v,{t)  = = + . 

The  integral  corresponding  to  the  first  term  is  estimated  by  ||ii'/i||//i/2  after  noting 
that,  for  the  relevant  values  of  s  and  <,  |5  +  t\/2S  <  1. 

The  integral,  corresponding  to  the  second  term,  can  be  estimated  by 


1 


/6'  J^  J^  \w,{t)\''dsdt  =  l/6\\w,\\l,^o^s), 


which,  in  turn,  can  be  estimated  appropriately  by  using  Lemma  6. 

The  third  diagonal  double  integral  is  estimated  in  exactly  the  same  way. 

We  next  estimate  the  off-diagonal  double  integrals.  We  note  that  for  0  <  <  <  (5, 
and  S  <  s  <  H  —  6, 

Ve{s)  -  Ve{t)  =  Wh{s)  -  -rWkit)  =  {Whis)  -  Wk{t))  + Wh{t). 

The  first  term  gives  an  integral  that  can  be  estimated  straightforwardly  in  terms  of 
|u^/i  1^1/2(0 //)•  What  remains  is  the  integral 


Jo      Js 


^^      *^     ds)\wHit)\Ut. 


6^{t-sy 


We  integrate  with  respect  to  s  and  find  that  the  inner  integral  is  bounded  by  1/6  and 
the  double  integral  by 

(23)  US\\w,\\l^o,sy 

The  estimates  of  the  other  integrals  can  be  carried  out  quite  similarly.  To  complete 
the  estimate  of  the  llvellui/Z/r^  zJ^■,  we  consider 


i 


( +  -7} )ds. 

s  n  —  s 


For  5  e  (0,(5),  and  5  G  {H  —  6,H),  we  obtain  contributions  that  Ccin  be  estimated  by 
the  expression  given  in  (23).  For  the  integral  over  ((5,  H  —  6),  we  use  Lemma  7. 

We  next  turn  to  the  space  associated  with  one  of  the  vertices  of  Q, .  We  now  use 
Sv{t)  =  (6  —  \i\)/6  to  complete  the  partition  of  unity,  i.e.  we  use  Ih{OyWh)  =  Ihi^v)- 
It  follows  from  Lemma  4  that  we  can  again  ignore  the  interpolation  operator  I^.  We 
need  to  estimate 


u 


-asat, 


s         \s-t\'^ 
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and 

J -8     6  —  s  s  +  6 

Considering  the  double  integral,  we  note  that 

{S  -  \s\)wKis)  -  {6  -  \t\)w,{t)  ^  w,is)-WH{t)      _  \s\  +  \t\    _  w,is)  +  w,(t)\s\-\t\ 
6{s-t)  s-t         ^  26      '  26  s-t   ' 

Since  |1  -  ^^^\  <  1  and  l^^fj^l  <  1,  for  relevant  values  of  s  and  t,  the  two  contri- 
butions to  the  integral  can  be  estimated  in  terms  of  |"'/>l^i/2(o,//)  ^^^  1/^II"^'»IIl2(o,5)' 
respectively.  Arguments,  quite  similar  to  those  given  above,  completes  the  proof  for 
the  case  of  two  dimensions. 

We  now  turn  to  problems  in  three  dimensions,  i.e.  the  case  where  the  interface  F 
is  two  dimensional.  In  addition  to  the  coarse  space,  we  use  three  types  of  local  sub- 
spaces  associated  with  faces,  Eind  neighborhoods  of  edges  and  vertices,  respectively. 
The  diameter  of  the  point  set  associated  with  a  vertex  subspace  is  on  the  order  of 
6.  Similarly,  the  edge  spaces  include  the  degrees  of  freedom  on  F  that  are  within  a 
distance  6  of  the  edge  in  question.  We  again  construct  a  partition  of  unity  associated 
with  these  sets.  As  before  these  functions  axe  continuous,  piecewise  linear  functions 
and  their  gradients  are  bounded  by  C/6.  The  proof  proceeds  as  in  the  case  of  two 
dimensions.  We  give  only  a  few  details.  We  use  ^e(^i)^e(<2)  to  construct  the  contri- 
bution to  the  decomposition  related  to  a  face.  Similarly,  we  use  6e{ti)dv{t2)  as  the 
part  of  the  partition  of  unity  associated  with  an  edge.  Using  our  formulas  for  9^  and 
9^,  we  then  show  that  the  partition  of  unity  is  completed  by  adding  functions  which 
differ  from  zero  only  in  small  neighborhoods  of  the  vertices.  The  estimates  necessary 
for  the  use  of  Lemma  1  and  the  completion  of  this  proof  are  then  carried  out  as  in 
the  two  dimensional  case.    D 
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